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<D ■ Abstract 

pq | We study 4d superconformal indices for a large class of M = 1 superconformal quiver 



gauge theories realized combinatorially as a bipartite graph or a set of "zig-zag paths" 



on a two-dimensional torus T 2 . An exchange of loops, which we call a "double Yang- 
lO ■ Baxter move" , gives the Seiberg duality of the gauge theory, and the invariance of the 

index under the duality is translated into the Yang-Baxter-type equation of a spin system 
CN ' defined on a "Z-invariant" lattice on T 2 . When we compactify the gauge theory to 3d, 

Higgs the theory and then compactify further to 2d, the superconformal index reduces 
to an integral of quantum/classical dilogarithm functions. The saddle point of this in- 
tegral unexpectedly reproduces the hyperbolic volume of a hyperbolic 3-manifold. The 
3- manifold is obtained by gluing hyperbolic ideal polyhedra in HI 3 , each of which could 
be thought of as a 3d lift of the faces of the 2d bipartite graph. The same quantity is 
also related with the thermodynamic limit of the BPS partition function, or equivalently 
the genus topological string partition function, on a toric Calabi-Yau manifold dual to 
quiver gauge theories. We also comment on brane realization of our theories. This paper 
is a companion to another paper summarizing the results pQ. 



1 Introduction and Summary 

Four- dimensional supersymmetric quiver gauge theories has been a useful playground to 
understand the physics of strongly coupled phenomena of gauge theories, in particular 
their IR fixed points. 

In this paper, we study a large class of 4d M = 1 quiver gauge theories described 
combinatorially by a configuration of oriented cycles (called zig-zag paths) on a two- 
dimensional torus, satisfying certain conditions analyzed below [2]. This combinatorial 
data, equivalently expressed as a bipartite graph (dimer) or a quiver diagram on T 2 , 
encodes the matter content (the quiver diagram) and the Lagrangian of our gauge theories 
[31 HI [5]. We here take the gauge group at each vertex of the quiver diagram to be U(N). 
The resulting gauge theory is believed to flow to a non-trivial interacting fixed point 
in the IR, can be engineered from a stack of N D3-branes probing the tip of the a 
toric Calabi-Yau manifold, and has been extensively studied in the context of AdS/CFT 
correspondence (see [U [7] and references therein) . 

There is an interesting subtlety in this story. The gauge theory corresponding to 
a given toric Calabi-Yau manifold is not unique, and several different gauge theories, 
related by a sequence of Seiberg dualities [8], correspond to the same geometry (this is 
sometimes called "toric duality"). In the language of zig-zag paths, this is translated into 
an ambiguity of the relative position of the loops, and Seiberg duality is translated into 
an exchange of the loops, which we call a "double Yang-Baxter move" . As the naming 
suggests, this is the standard Yang-Baxter move repeated twice, and strongly suggest the 
integrable structure behind the theory. 

Given a 4d supersymmetric gauge theory, we could extract concrete quantitative 
statements of the theory by computing its 4d superconformal index / (9j [TO]. This is a 
twisted partition function on S 3 x S 1 , where the chemical potentials are turned on along 
the S l direction. The index can be computed in the free field limit, and is written as a 
matrix integral. One of the main results of this paper is that this matrix model could 
be regarded as the partition function of a spin system defined from zig-zag paths on T 2 , 
where each spin has N — 1 continuous values in S 1 . This is summarized in the relation 

■L 4d quiver — ^ spin system on T 2 ■ \ *-•*- ) 

This is the manifestation of the integrable structure mentioned above; the Seiberg dual- 
ity is now translated into the statement that the resulting partition function is invariant 
under the double Yang-Baxter move, ensuring the integrability of the model. Interest- 
ingly, (modulo some important differences mentioned below) the resulting spin system is 



essentially the same as the spin system studied in [TT] for SU(2) gauge groups, and more 
recently in [12] for SU(N) gauge groups. 

Let us next study the reduction of our 4d M = 1 theories along S l . The resulting 
theory has 3d M = 2 supersymmetry, and flows in the IRg to a non-trivial fixed point. 
As the S 1 shrinks all the KK modes decouple, and the 4d superconformal index should 
reduce to a partition function on S 3 . Indeed, it has been shown that the 4d index in 
this limit reduces (after suitably regularizing divergences) to a 3d partition function on 
ellipsoid S% (defined in (j4.ip ). which could be again written as a matrix integral after 
localization computation: 
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This limit is also natural in the context of integrable models; the Yang-Baxter equation, 
being an equality, should hold even after taking the limitp. After taking one more limit 
explained in the text (Higgsing to the Abelian gauge group), the solution of the star- 
triangle relation studied in [TT] reduces to another solution discovered by Faddeev and 
Volkov [13], dH [15] , clarifying the integrable structure behind 3d M = 2 theories. 

We could also consider further dimensional reduction to 2d. This is simply the b — > 
limit of the ellipsoid partition function, and taking the leading contribution we have 



^3d on Sl ► ^2d on R 2 = / d <7 exp 
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where a is the scalar component (s) of the twisted superfield (defined from the derivative 
of the vector superfield) and takes values in the Cartan of the gauge groupo The poten- 
tial H ; 2d( cr ) represents the effective twisted superpotential obtained by integrating out 
matters from the theory. 

The surprising observation, based on the works [TBI [T7], is that this twisted effective 
superpotential is identified with the hyperbolic volume of a certain 3-manifold M, in the 
case that N = 2. 

The 3-manifold M is determined from the bipartite graph on T 2 which in turn is 
determined from zig-zag paths, and could be thought of as a 2d graph with an "extra 
dimension" added. The 3-manifold M is defined as the union of ideal hyperbolic polyhe- 
dra in HI 3 , and the projection of the polyhedra onto the boundary of H 3 gives the faces 
of the 2d bipartite graph. 



1 As we will see there are some subtleties associated with the regularization of divergences in the limit. 
2 In general a is a vector, but we here do not show this fact explicitly for notational simplicity. 



The twisted superfield scalars a, in this description, is identified with the geometric 
modulus of the 3-manifold M; the dihedral angles of M are determined from the radii 
of circles on T 2 , whose logarithms coincide with a. The values of a are determined from 
the gluing conditions of the 3-manifold M 

fdVo\[M}(a) 
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There is a counterpart of this equation on the gauge theory side; the value of a at the 
vacuum is determined from the equation 

'0W 2d 



cxp 



da 
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We find that the two conditions (jl.4p . (jl.5p coincide. In other words, the vacua of the 
2d M = (2, 2) theory is captured by the gluing condition of the 3-manifold! This is the 
second main result of our paper. 

One quick supporting evidence for the correspondence between (jl.4p and ( II. 5p is that 
the twisted superpotential W2d is expressed as a sum of the Lobachevsky functions (or 
classical dilogarithm functions), and the same function is known to appear in the formula 
for the volume of hyperbolic tetrahedra; M is simply the sum of these tetrahedra. Of 
course, the appearance of the dilogarithm function applies to any 3d M = 2 theories 
dimensionally reduced on S 1 , whereas our correspondence should hold only for a specific 
class of 3d gauge theories. 

In Table [I] we summarize our correspondence between 2d M = (2, 2) theory and the 
geometry of the 3-manifold M. The data on both sides come from the zig-zag paths on 
T 2 , and therefore from a toric Calabi-Yau 3-fold or from the brane configuration for our 
gauge theories. 



Table 1: Dictionary between the 2d M = (2, 2) theories and the 3-manifold. 



2d jV = (2,2) gauge theory 


3-manifold 


twisted superpotential W2d(c) 


Vo\[M](a) 


scalar in twisted superfield a 


modulus a of M 


matter contributing Li 2 


tetrahedron contributing Li 2 


vacuum equation exp (^r) = 1 


gluing condition exp (^r) = 1 



The correspondence to this point refers only to the 2d gauge theory. However, it is 
natural to ask if similar correspondence persists for the 3d/4d gauge theories we started 



with. As for the 3d gauge theory, the natural guess is to propose 



J 3d SL(2) Chern-Simons on M j 
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where the right hand side is the (holomorphic) partition function of the SL(2) Chern- 
Simons theory on M and the parameter b is identified with the inverse square root of the 
level t of the Chern-Simons theory: b 2 ~ l/(t + 2). This is consistent with our previous 
correspondence since the classical limit of the SL(2) Chern-Simons theory reproduces the 
volume (and the Chern-Simons invariant) of the 3-manifold. There is also generalization 
of (II. 6p to N > 2, where the right hand is replaced by the partition function of SL(N) 
Chern-Simons theory on M. 

The relation fll.6|) . relating 3d N = 2 quiver gauge theories and the 3d SL(2) Chern- 
Simons theory, is highly reminiscent of the recently found connection [T8l [191 l20l |2"T] |22| 
[231 121] between 3d M = 2 theories and 3d SL(2) Chern-Simons theories (see also |25j). 
There one of the crucial underlying data is the Riemann surface and its Teichmiiller 
space, whereas here we have a dimer as the crucial ingredient. We expect that this 
similarity could be explained from the equivalence of underlying mathematical structures 
(for example, cluster algebras), perhaps along the lines of [26] (c/. [2TJ, [23], [28]). It would 
be interesting to understand the precise relation between the two. 

Let us also comment on the relation between (11.11) and ( 11.61) . It is natural to interpret 
both statements from compactification of a 6d theory (see the discussion in section 12.31 
and 14.51) . If this is true, then the two statements are related by a dimensional reduction 
on the one hand, and by a dimensional oxidation on the other side, thus exemplifying 
the statement (c/. (29] section 5]) 

dimensional reduction = dimensional oxidation 

in the AGT[30]-type correspondence. 



Finally, we point out connection of our results to topological string theory and the 
BPS state counting. 

Under an assumption about the bipartite graph (isoradiality condition in section |2~2]) . 
we show that the critical value of the hyperbolic volume of our 3-manifold M could be 
written as a sum of Lobachevsky functions (see ( 14.38P ). Interestingly, exactly the same 
expression arises as a Legendre transform of the thermodynamic limit of the partition 
function of the dimer model. This dimer model has been studied in the context of BPS 
state counting of type IIA string theory on a Calabi-Yau 3-manifold X&, which in turn is 



known to be equivalent to the topological string partition function (modulo wall crossing 
phenomena). In this context the thermodynamic limit is the semiclassical limit g top — > 0, 
where g top is the topological string coupling constant, and the leading contribution is 
precisely the prepotential J-to P ,o [3T] . Combining these observations, we have 

Ftopfi is an integral of the Legendre transformation of Vol[Mo] . (1-7) 

This paper is companion to [1], which announces basic results. 

This paper is organized as follows (see Figured] for the logical structure of this paper). 
After a summary of 4d M = 1 quiver gauge theories described from zig-zag paths (section 
|2]), we compute its superconformal index and comment on the reformulation as a spin 
system (section [3]). We then reduce the theory down to 3d and 2d, and study the 
connection with 3-manifolds (section H]). Section [5] explains the relation with topological 
string theory and the statistical mechanical model of BPS state melting. We conclude 
with some future problems (section [6]). Appendix contains a summary of the special 
functions used in the main text, and an explicit computation of thermodynamic limit of 
the dimer partition function. 
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Figure 1: Logical structure of this paper. Clearly it is impossible to list all the connections between all the ingredients 
mentioned here. The main claims of this paper are the two equalities represented in the center of this figure. 



2 Quivers from Zig-Zag Paths 

In this section we briefly summarize the construction of 4d M = 1 superconformal quiver 
theories from "zig-zag paths" on T 2 (2j [32] (see (33j [26] for mathematical formulation). 
See also the reviews [6j [7] for more details on dimer model techniques. 

2.1 Zig-Zag Paths 

For the clarify of the presentation let us first explain the combinatorial properties of the 
zig-zag paths, which is actually rather elementary. The physical context will be explained 
shortly. 

Let us start with a convex polygon A in Z 2 . Geometrically this is the toric diagram 
for a Calabi-Yau 3-manifold X&, i.e., the cone A x {1} G IR 3 specifies the top-dimensional 
cone of the fan. 

As a toric diagram there are SL(3, Z) ambiguities in the choice of A. For example, 
Z 2 translation of A keeps the geometry. In the following we use the same symbol A for 
the equivalence class of A under this identification. 

One way to specify A is to write down the set of primitive normals of the polygon. 
Let us denote them by (r^Si) ^ (0,0) with i = l,...,d. We choose the label i such 
that the direction of the vector p^ rotates in the counterclockwise manner as we increase 
i. The integer d, which is the number of lattice points in the boundary of A, is fixed 
throughout this paper. Note that in general the same vector could appear multiple times 
in this list. This happens when an edge at the boundary of A contains more than two 
lattice points. 



By definition we have 



$> = 5> = 0. (2-1) 



Let us now consider zig-zag paths. The zig-zag paths are a set of closed oriented 
cycles p\, . . . ,pa on a two-dimensional torus T 2 , whose homologies cycles are determined 

by (ri,Si): 

[Pi] — r A a ] + Si\P\ £ Hi(T 2 ,Z) fixed and non-trivial , (2.2) 

where [a], [/3] are the basis of Hi(T 2 , Z), for example a and /3-cycles of the torus. There 
is SL(2, Z) ambiguity in the choice of [a], [ft], which could be absorbed into the SL(2, Z) 
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ambiguity in the choice of A. The origin of the terminology "zig-zag path" will be 
clarified shortly when we study dimer models. 
We assume the following three conditions. 

• genericity. First, we assume that no three paths intersect at a single point (Figure 
[2]). This is satisfied for generic choice of paths. 

• admissibility. Second, we impose the admissibility condition (this terminology 
comes from [33]). To explain this, let us note that the paths divide the torus into 
a union of convex polygons. We color the convex polygon by black (white) if all 
the paths around the polygon has counterclockwise (clockwise) orientation around 
the polygon; otherwise the face is kept uncolored. The paths are called admissible 
if every edge bounds a colored polygon (Figure [3]). 





Figure 2: Genericity condition, stating that no three zig-zag paths intersects at a single 
point. The left figure is allowed, whereas the right is not. 

These two conditions are sufficient for the 4d quiver gauge theory. We moreover 
impose one simplifying assumption 

• minimality. Minimality (this terminology comes from [26]) forbids the two possi- 
bilities shown in Figure HI 

We will come back to the physical significance of the admissibility condition in section 



12. 3[ but for the moment let us first analyze its combinatorial implication by defining 
graphs on T 2 . 

Given a set of zig-zag paths we can define a natural bipartite graph Q* and its dual 
Q, both realized on T 2 (the symbol is chosen for later convenience). In the literature 
Q* is often called a brane tiling [31 IU [5], and Q a periodic quiver. The vertices of Q* 
are given by colored faces, and the edges by the intersection points between them. The 
orientation of zig-zag paths gives a natural orientation to the edges, and ensures that 
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admissible 



non-admissible 





Figure 3: Admissible (left) and non-admissible (right) configuration of zig-zag paths. 
Rather than coloring the faces by black and white, we have represented the coloring by 
placing black and white dots inside (it is hard to represent the white color on a white 
paper!). We see from this example that moving a zig-zag path across an intersection 
point breaks the admissibility condition. 





Figure 4: Minimality condition forbids two types of intersections of zig-zag paths. The 
graph superimposed on it is the bipartite graph Q*. 



the resulting graph is bipartite, i.e, vertices are colored either black or white and edges 
connect vertices of different colors. 

The graph Q* is dual of Q. This means the vertices are placed on the uncolored 
faces, and the edges at each intersection point between them. The orientation of the 
zig-zag paths determine the orientation of the edges. In the following we denote the set 
of vertices, edges and faces of Q by V, E and F. Since we have black and white colors 
we have a decomposition into black and white faces: F = B U W. For an edge e e E we 
denote the source (and the target) by s(e) (t(e)). Since Q is drawn on T 2 we have 



V 



E 



. 



(2.3) 



We also denote by V*, E*, F* the set of vertices/edges/faces of the graph Q* . By definition 
we have 

V* = F, E* = E, F* = V . 

We have defined Q, Q* from zig-zag paths, but we can go in the other direction. Given 
a bipartite graph we define a zig-zag paths to be graph on Q* which turns maximally 
right (left) at black (white) vertex. Because the graph is finite, we always come back to 
the same vertex after several steps and hence this defines a set of closed loops. The name 
zig-zag path originates from the zig-zag shape of the path_| (see Figure [5]). 




Figure 5: A zig-zag path on the bipartite graph Q* (dotted path) is identified with the 
zig-zag path defined previously (undotted arrow). 



3 This is also called a rhombus loop or a train track in the literature. The word "rhombus" refers to 
a quadrilateral in Figure [TU1 which becomes a rhombus for isoradial circle patterns. 
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We can verify that this gives an inverse to our previous construction, i.e., from the 
bipartite graph we recover the zig-zag paths we started with. However, it is important 
to notice that this correspondence is not one-to-one; we can start with an admissible 
configuration to obtain another configuration by locally applying the two basic moves 
shown in Figure |6j These two moves are called fundamental moves in this paper and 
complete in the sense that any two minimal bipartite graphs corresponding to the same 
A are related by a sequence of these two moves [26l Theorem 2.5]. The same theorem 
guarantees the existence of minimal admissible configuration of zig-zag paths. Note that 
the number of nodes of the quiver is preserved in the fundamental moves. 



Move II 



Move I 





Figure 6: Two basic moves (fundamental moves) preserving the admissibility condition. 
The second move is called a double Yang-Baxter move. 

It should be kept in mind that the choice of the fundamental moves is not unique. For 
example, we could replace the second move in Figure [7] by a different move, for example 
the ones shown in Figure EJ We can easily check that the these moves, in combination 
with move I, generate the same set of moves. We call move II' a Seiberg move (it is exactly 
Seiberg duality, as we will see shortly) and move II" a spider move. For later reference, 
we also list basic moves for zig-zag paths without admissibility condition imposed (Figure 
E|) • These are reminiscent of the Reidemeister moves of knot theory. 
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Move II' 





Move II" 





Figure 7: We can replace the fundamental move II in Figure |6] by either of the two moves 
shown here. 







Figure 8: Moves for zig-zag paths without admissibility condition. Here we do not show 
the orientation of the zig-zag paths, and all orientations are allowed. 
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2.2 Quiver Gauge Theories 

Let us next define 4d M = 1 quiver gauge theories from the combinatorial data of the 
previous subsection. This is simply a quiver gauge theory determined from Q. In other 
words, we have a gauge group SU(N) V for each vertex v G V, and a bifundamental chiral 
multiplet X e for each e E E. The total gauge group is given by 

G = J] SC/(iV)„ . (2.4) 

vev 

Here and in the following we will specialize to the case where the ranks of the gauge 
groups are all equal to N. Most of the analysis of the next section on 4d superconformal 
index generalizes straightforwardly to more general cases where the ranks are position 
dependent. However, the analysis of the Seiberg duality in 4d, and also of the reduction 
to 3d, requires some change. 

We also determine the superpotential to be 

W = Y,^(U X )-J2^(U X ) > ( 2 - 5 ) 

b£B \eeb / weW \e€«) / 

where the product inside the trace is taken in the counterclockwise (clockwise) manner 
for B and W. 

The quiver gauge theory constructed in this way contains examples with enhanced 
supersymmetries, for example the M = 4 theory or M = 2 theories corresponding to 
the Al type singularities [34]. However, generically the theory is chiral, and has M = 1 
supersymmetry. The typical example in the literature is the theory dual to the conifold 
[35]. 

Although our theory is in general chiral, there is no chiral anomaly. This is because 
by construction the number of incoming and outgoing arrows for each vertex v G V are 
the same. 

The quiver gauge theories constructed in this way is the world- volume effective theory 
on the D3-branes probing the toric Calabi-Yau manifold X&, and flows in the IR to a 
non-trivial IR fixed point. 

The basic moves in Figure |6] is an operation on quiver gauge theories which keeps the 
IR fixed point intact [I]. The first move is to remove two fields X, Y with superpotential 
Ty(XY) from the theory — because this superpotential term represents a mass term, we 
can simply integrate out a massive field. The second move (move II') corresponds to a 
Seiberg duality (see [7], section 4.7 for detailed exposition), or equivalently the mutation 
of the quiver. 
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Note that in our setup Seiberg duality can be taken only for those nodes which has 
four arrows (two incoming and two outgoing). This is because Nf is always a multiplet 
of N, and the only value of Nf in the conformal window is Nf = 2AO- 

Finally, let us comment on a further simplifying assumption on the bipartite graph. 
A bipartite graph is called isoradial if all the vertices can be placed on circles of equal 
radius. In terms of zig-zag paths this is satisfied if and only if (1) each zig-zag path 
is a simple closed curve and (2) the lift of any pair of zig-zag paths to the universal 
cover intersect at most once [36] • Colloquially this means that all the zig-zag paths are 
"straight enough" such that the zig-zag paths have minimal intersection numbers. This 
means that the number of edge of the quiver is minimal, i.e., we have 



\e\ = J2 \faPi)\ = £ 



r i S j 



r jSi\ 



(2.6) 



Kj 



i<j 



We will see that this condition leads to enormous simplifications of part of the upcoming 
analysis; for example, in section 14.31 we will comment on this condition in the context of 
2d gauge theories. However, it should be kept in mind that isoradiality is not a necessary 
condition, and for example stronger than the consistency conditions in [32} [371 EH] • 

2.3 Brane Configuration 

The quiver gauge theories in the previous subsection could be realized by D5/NS5 brane 
configurations [39l HOI [7]. Let us briefly summarize this, since this clarifies the origin of 
the apparently ad hoc assumptions in the previous subsection. 
The relevant brane configuration is 



(2.7) 



Here N D5-branes wraps the T 2 along the 57 directions, whereas a single NS5-brane 
wraps a holomorphic cycle E in 4567 direction. If we write x = e Xi+lX5 ,y = e X6+ * X7 , then 
the holomorphic curve is given by 








1 


2 


3 


4 


5 


6 


7 


8 


9 


7VD5 


- 


- 


- 


- 














NS5 


- 


- 


- 


- 




5 


i 









x\2 



{P(x,y) = 0}c(C x ) 



(2.8) 



It is important that not all the possible mutations of the quiver physically make sense. If we mutate 
on the n-valent vertex with n > 6, the result is in general a graph which could not be realized on torus. 
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where P the so-called Newton polynomial of A, defined by 

P(x,y)= Y, C ^V • (2-9) 

(»,i)eA 

with generic coefficients cy. Here (z, j) G A means the lattice point (i,j) is a lattice 
point of A (including the boundary). The curve E is a Riemann surface of genus g and 
puncture d, where g is the number of internal lattice points inside A (recall that d also 
denotes the number of zig-zag paths in section [2TTT) . This curve is part of the mirror X& 
of the toric Calabi-Yau 3-fold X& given by [UJ 02] 



X A :uv + P(x,y) = , (2.10) 

with u, v G C. We will encounter P(x, y) again as a spectral curve of the dimer model. 
The identification of the NS5-brane curve and the mirror curve could be explained by a 
T-duality, where our D5/NS5 system is mapped to a configuration of D6-branes wrapping 
Lagrangian 3- cycles inside X& [4"3] . 

This brane configuration can be looked at from two different ways: D5-brane view- 
point and the NS5-brane viewpoint. In the former, we have a T 2 , which is the torus we 
had previously, and the zig-zag paths represent the intersection cycles of an NS5-brane 
and the D5-brane; and the colored faces represents the projection of the shape of the NS5- 
brane. In this language, the colored regions correspond to the projection of the mirror 
curve into the 57 directions (called coamoeba/alga [12]), and the black/ white color rep- 
resents the orientation of the curve when projected onto T 2 . In other words, black/white 
regions represents (N, 1)/(N, — l)-branes, and uncolored regions (N, 0)-branes. This ex- 
plains why U(N) gauge group lives in uncolored regions. The bifundamental fields origi- 
nate from the massless strings between the U(N) gauge groups, namely the intersection 
points of the uncolored regions. The admissibility condition simply says we do not have 
(N, fc)-branes with \k\ > 2, in which case no Lagrangian descriptions are known. 

We could also take the NS5-brane viewpoint. Then we have a Riemann surface E, 
on which we have a set of 1-cycles representing the intersection with D5-branes. Our 
Riemann surface E is reconstructed as the surface whose boundaries are the cycles of zig- 
zag paths jl3] . This is parallel to the construction of the Seifert surface in knot theory, 
and Seiberg duality in the original graph is translated into the half Dehn twist of E [26] . 
We will encounter this Riemann surface E again in section 14.51 
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2.4 R-charges 

We are going to consider the IR fixed point of our quiver gauge theories. Due to the 
strong coupling effects the bifundamental chiralmultiplets could have large anomalous 
dimensions in the IR. This is determined from the IR U(1)r R-symmetry in the su- 
perconformal algebra. In general, it is a rather non-trivial problem to identify the IR 
superconformal R-symmetry U(1)r, because UV R-symmetry mixes with the global UV 
U(l) symmetries. Here we will comment on one particular useful parametrization of the 
UV global symmetry or equivalently IR R-symmetry (see [21 EH1 Pf0) 5 l ). 

We are going to define d global symmetries, with 1 relation among them. For each 
zig-zag path p t , let us define the charge of the bifundamental field X e at an edge e G E 
by 

Qi[X e ] = <pi,e), (2.11) 

where the braket here refers to the (signed) intersection number of the two paths Pi and 
e. To see that this is in fact a global symmetry, recall that a term in the superpotential 
is represented by a closed loop around a black/white vertex of the bipartite graph (see 
(12. 5p ). and has intersection number with a closed loop p- h . Because two zig-zag paths 
pass through e with an opposite orientation, we find that the diagonal subgroup of these 
d global symmetries is trivial (see 12.11) 



£)Qi[X e ]=0. (2.12) 



Hence we find d — 1 global symmetries. These symmetries are anomaly free [H] , because 

J2 Qi[X e ] = for alii. (2.13) 

e: around v 

This holds because a zig-zag path incoming to a vertex necessarily goes out of the vertex, 
with an opposite orientation and hence with a field with an opposite flavor charge. 

It is known that generically this parametrization exhausts all the possible global 
symmetries^!. There is a corresponding statement in the AdS dual; two out of d — 1 
correspond to isometries of the Sasaki-Einstein manifold and called mesonic symmetries, 
whereas the remaining d — 3 symmetries are associated with the 3-cycle of the Sasaki- 
Einstein manifold and called baryonic symmetries. 



5 The paper [2] discuss the case of isoradial bipartite graphs, however our parametrization applies to 
more general bipartite graphs. This will be crucial when we discuss Seiberg duality in section ["3.21 
6 In some speical cases there could be an enhancement to non-Abelian global symmetries. 
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We can also describe R-symmetries. The choice of UV R-symmetry is not unique, 
since we can consider a mixing with the d — 1 global symmetries mentioned above (the 
IR R-symmetry inside superconformal algebra is determined by a- maximization (35]). 

For our purpose, a particularly useful parametrization is given as follows. Let us 
choose a set of d parameters #, for each zig-zag path p^, and let us assume that they are 
defined modulo 2tt ($i ~ 0» + 2ir), and that < 9 i+i — 9 t < tx for all i (9 d+1 := 9 ± ). We 
can regard 9i as the slope of pi. Then for each bifundamental X e for e G E its R-charge 
is simply defined to be the relative slopes of the two zig-zag paths which goes through 
the edge. More formally, we define 

R e = R[X e ] := -sign(p L{eh p R{e) )[9 L{e) - 9 R{e) ] , (2.14) 

where [x] denotes a real number in [0, 2ir] and equivalent to x modulo 2ir. We also used 
the notation that a chiralmultiplet corresponding to an edge e has a flavor charge +1 
for the L(e)-th zig-zag path and —1 for the R(e)-th path, where Pi{e) and PR( e ) are two 
zig-zag paths passing through e from opposite sides (see Figure [9]). From the definition 
we have R P > 0. 




Figure 9: An bifundamental at an edge e has flavor charge +1 for L(e)-th flavor charge 
and —1 for i?(e)-th flavor charge, where PL(e),PR(e) are two zig-zag paths as in this Figure. 

We have the following two conditions on IR R-symmetries. 

• First, the /3-function for Yukawa couplings vanish. This is the same as the re- 
quirement that the R-charge of the superpotential, and therefore any term in the 
superpotential, is normalized to be 2. This means 

Y,Re = 2. (2.15) 

eS-F 

• Second, the /3-functions for the gauge coupling vanish. From the NSVZ /3-function, 
which in our case could be written as 

d 1 N 



dlogiigl 



glN/Sir* 



5 Ed 



7e 



e£v 



(2.16) 
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where the anomalous dimension 7 e is related to the R-charge by 7 e = 3R e — 2. 
From this condition, we have 

X>- Re) = 2. (2.17) 

e£V 

This can be written more symmetrically 

Yl R *e = 2 > ( 2 - 18 ) 

eGF* 

where we defined R* := 1 — R e . 

These two conditions follow from the definition (12. 14ft and the fact that the sum of 
exterior angles of a polygon is 2ir. 

Naively the dimension of the solution space to (1 2 . 1 5 j) . (12.17]) is zero because we have 
\E\ parameters R e and \V\ + \F\ constraints (see (J2.3P ). However, in supersymmetric 
theories not all the constraints are independent and it has been demonstrated that the 
solution has d — 1 parameters [ID]. Again, it is straightforward to see that the two 
conditions (I2.15p . (12. 18[) are preserved under the deformation with the d — 1 parameters. 
Of course, this is the same as the number of global symmetries treated above. 

3 Superconformal Index As Spin System 

In this section we first define the superconformal index for 4d M = 1 superconformal 
gauge theories. We then show that the 4d index for the spin system defined in previous 
section is equivalent to the partition function of a spin system on T 2 . 

3.1 Superconformal Index 

The Definition 

Let us consider 4d M = 1 superconformal theory on S 1 x S 3 . This theory has 
supercharges Q a , Qa and S a ,Sa, where a and a denotes the spins SU (2)i and S , f/(2) 2 
of the isometry of S 3 : Spin(4) = SU(2)i x SU(2) 2 - To define an index, we need to pick 
up a particular supercharge. There are four supercharges, but Q + and <2_ (Q_j_ and Qj_) 
define the same index due to the SU(l)i (5 , 6 r (2) 2 ) symmetry, and hence we have two 
possibilities. When we choose Q = <2_, the superconformal index is defined as an index 
defined from Q, with insertions of operators commuting with Q: 

I L (t,y; z) = Tr [(-1)^+^^-^^] , (3.1) 



where the index is taken over the Hilbert space on S 3 . F is the fermion number, £ is the 
energy, j x and j 2 are the spins of 577(2) i and S77(2) 2 , respectively, and u T ■=\\ i u i % 
where Ti (ui) is the charge (chemical potential) under the i-th. flavor symmetry^]. Note 
that due to the translation symmetry we could flip the sign of ji, and we have 

I L (t J y; U ) = I L (t,y~ 1 ;u) . (3.2) 

This index, often called the left-handed index (hence the symbol L in (13.11) ) is independent 
of the value of 7 thanks to the standard index argument, and can be computed in the 

limit 7—7-00 

I\t, y- u) = Tr [(-1) V^+^tf*] , (3-3) 

where the trace is now taken over all states satisfying 

{Q,&} = £-2j 2 + ^r = 0, (3.4) 

where r is the U(l) R-charge. Similarly, if we choose Q = Q^, we can define the right- 
handed index 

I R (t, y; u) = Tr [{-1) F t^ E+j ^y 2J2 u F ] , (3.5) 

where the trace is taken over all states satisfying 

{Q,Q J } = E-2 Jl -K = 0. (3.6) 

This is the same as the left-handed index, except that the orientation of the arrows of 
the quiver diagrams are reversed. Since this can be taken into account by a change of 
convention, in the following we will concentrate on the left-handed index and denote the 
corresponding index simply by I. 

For our later purposes it is useful to reparametrize the chemical potentials as 

p = t 3 y, q = tV 1 • (3.7) 

In this notation, the superconformal index reads 

I L (p, q; u) = Tr [(-l^p^+^g^-^e-^'^l . (3.8) 

One subtlety we have is that our 4d quiver theory is defined in UV, and is conformal 
only in the IR, where the theory is strongly coupled. Because the index is independent 



7 This is often denoted by z.i in the literature. 
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of the continuous parameter, it is independent of the dimensionless parameter obtained 
by multiplying the energy scale by the radius of S* 3 , and we could compute the index in 
UV, except that we have to take into account the mixing of UV R-symmetry with global 
symmetry. In our analysis this is taken into account when we have included chemical 
potentials for global symmetries. In practice this means that the effect of anomalous 
dimensional could be taken into effect by shifting the global symmetry chemical potentials 
by powers of t. 

Integral Expression 

Because an index is invariant under the continuous deformation of parameters of the 
Lagrangian, the superconformal index can be computed by taking the free field limit. 
Alternatively, we could apply localization techniques. In either way, the result is written 
as an integral over the Cartan of the Cartan H of the gauge group, which in our case is 
given by 

H = (l^l)^ 1 ) lVl CG = SU(N)W . (3.9) 

We parametrize an element of the Cartan of SU(N) V (at vertex v G V) by z v = 
(z Vt i, . . . ,z Vt x) satisfying Yli=i z v„i — 1- We also write z V j = e lCTv >% where a v>i is peri- 
odic with period 2n. Physically these parameters represents the Polyakov loop along the 
thermal direction S 1 . 

The index is expressed in a plethystics form (this follows from group theory, see [ID] ) 



I(p, q,u) = Yl 



dfi v Y[(z vi - z vJ )(z v> l - z„j) 



exp (^2i{p n ,q n ,u n ;z n )) (3.10) 



Here the integration measure contains dfx v defined by 

AT-l , 

az 



,n=0 



*. = ^II£r. ( 3 -"= 



iV! J-J- 27TZZ 



as well as the Vandermonde determinant, and the integral over the contour \z v ^\ = 1 
(or a Vt i runs from to 2n). The "single-letter index" i(p,q,u;z) is given as a sum over 
contributions from vector and chiralmultiplets 

l(p, q, U]Z) =^ KectiP, Q'iZ) + ^2 ^chiralO; Q, u \ z ) , (3-12) 
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where 

1 — pq 



%ect(p,Q',z) 



1 - 



l-p)(l-q) 

1 

l-p)(l-q) 



Xadj(^), (3.13) 



1 [Re 

^chiralfe <?, "I *) = n _ V1 _ x (pg)"^y«L(e)«fl( e )Xbiftind(^(e), ^t( e ); 



- (P?) 1 2 7(Z/WL(e)M if ,( 1 e) )Xbifund(2; s( e ) ,^( e ) 1 )J , (3-14) 

and Xadj and Xbifund are the characters for the adjoint and bifundamental representations, 
respectively 

N 

Xadj(z) = ^ i Z i Z j l ) ~ h Xbiiund(z,z') = ^ ^j 1 ' ( 3 " 15 ) 

l<i,j<N i,j=l 

Note that vector multiplets are not charged under flavor symmetries. 
For our purposes, it is useful to rewrite the index in a different form 

I(P, q,V>)= j[[ d ^ v II J vcct(P 5 q\z)\ f XI 7 <*iral(P> 9, u; z) 



where 



C*(p, ?; *) = *fe qf' 1 II f? — ^ — v ' ( 3 - 17 ) 



and 



^hiralfe 95 *) = II r ((PS0 ^ «i(e)Wfl(e)^(e),fc^(e),i;P» «) " ( 3 ' 18 ) 

l<fc,Z<iV 

Here r(z;p, g) is the elliptic gamma function defined in (lA.lj) and we defined (see (IA.5I) ) 



K(p,q) := {p;p)oo{q;q)oc ■ (3.19) 

The equivalence of the two expressions (I3.10p and (13. 16[) can be verified by the equa- 
tions of the form 



ex (v 1 Vn 1 

\n=l / m=0 



(3.20) 



and the equalities in Appendix, for example ( 1A.8J) . 

We here obtained the expression of the index (13. 16[) by rewriting (I3.10p . However 
it should be emphasized that ( 13. 16[) arises directly in the localization derivations of the 
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index (see for example [16] and [29j Appendix]). The factors I^ ect (/matter) represents 
the 1-loop determinants for vector multiplet at vertex v (chiralmultiplet at edge e), and 
the infinite product comes from the spherical harmonics expansion on S 3 . In the free 
field computation, we have the S 3 Laplacian from the bosons, whose determinant could 
be written as products over the quantum numbers of the spherical harmonics, and the 
same applies to fermions. There are cancellations between bosons and fermions, and 
the unpaired bosonic (fermionic) modes appear in the denominator (numerator) of the 
product flA.ljl . 

3.2 Invariance under Fundamental Moves 

In this section we prove an invariance the superconformal index under the fundamental 
moves. The general argument that the 4d index depends only at the IR fixed point 
guarantees this invariance, but we can also check this explicitly. 

Generalities on Gluing 

It is important to note that the two fundamental moves in Figure[6]are local operations 
on the quiver diagram. We therefore expect that the invariance of the index should reduce 
to the invariance of the index defined for the sub diagram. 

To formalize this idea we need to invoke the concept of "gluing" in gauge theories. 
Let us first explain this in a rather general situational Consider two 4d N = 1 theories 
T\ and T 2 which has global symmetries G\ and G 2 . Suppose moreover that these flow in 
the IR to non-trivial fixed points. We can then compute the superconformal indices for 
each of these theories, with the chemical potentials for global symmetries included. 

Now from the two global symmetries G\ and G 2 we choose a common subgroup 
H (Hi C Gi, H 2 C G 2 , H ~ Hi ~ H 2 ) and gauge H, i.e., the diagonal subgroup of 
Hi x H 2 . The resulting theory has global symmetry Gi\Hi x G 2 \H 2 , where G\H denotes 
the commutant of H inside G, see Figure [TUl 

This gluing operation has a counterpart at the level of the index. Let us denote the 
indices of Ti and T 2 by Ii(p,q;ui,w) and I 2 (p, <?; u 2 , w), where Ui and u 2 denotes the 
chemical potentials for the Cartan of the flavor symmetries Gi\Hi and G 2 \H 2 , respec- 
tively, and w for those for H. We then havqj 

—— / Ii(p,q;ui,w)I vect (w)I 2 (p,q;w,U2) . (3.21) 



8 This discussion obviously generalizes to quantum field theories in other dimensions, for example to 
ellipsoid partition function in 3d studied in section T4. II 

9 This is the gluing of "generalized index" , see [27] for similar analysis in 3d. 
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gauge j f ungauge 




Figure 10: We can glue two theories T\ and T^ by gauging a common global symmetry 
H. 

Note that after gluing the global symmetry is promoted to dynamical degrees of freedom 
to be integrated over, and hence we have included the 1-loop determinant for w. 

Let us apply this to our case, where one of our moves replaces a quiver diagram 
with another. Since the fundamental moves are local operations on the graph, we could 
decompose the quiver diagram T into two parts along edges of T such that Ti becomes 
T2 after the move and T\Ti = T\T2 kept intact. Then the rule in ( 13.211) immediately 



gives 



and 



J r = / dfx(z) I vcct (p, q) I Tl (p, qi z) I r/Tl (p, q; z) , 
J r , = / dfi(z) IvectO, q) ir a (p, q; z) I r /r 9 (P. T, z) . 



(3.22) 



(3.23) 



The equality of these two quantities follow if we could show the local invariance of the 
index 



-fri - h 2 ■ 



(3.24) 



Because any duality is generated by combination of the two moves (section I27L]) . all we 
need to do is to check the invariance under the two fundamental moves. 

Integrating Out 

Let us first study move I. After the move we have two extra bifundamental fields 
X\ and X 2 between two gauge group SU(N) V and SU(N) V >, which has a superpotential 
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term Tr(XiX2). This means that the sum of R-charges is two and their global symmetries 
charges have opposite signs, and the contribution from the two is 

T {(pq^zvzjuiufwq) T ((pqf-^^z^z^ufu^p^) . (3.25) 



But this is trivial due to (1A.2J) . 

Seiberg Duality 

Let us first analyze move II', which is equivalent to the Seiberg duality The invariance 
of the index under Seiberg duality has been verified in jlHj. However, an extra analysis is 
required here because the invariance of the index holds only for particular assignment of 
R-charges to fields, and we need to check that the assignment of the R-charge in section 
12.41 satisfies this condition. 

As we demonstrated already we can concentrate on the part of the bipartite graph 
which changes under the duality. Let us label the fields of the electric theory by Xi and 
the magnetic theory by Y iy Z t (i = 1, . . . , 4) (see Figure [11]). Following the rule (j2.14p 
and using the quantity 

Rij := -83gn(j>i,pj) [9i - 6j] , (3.26) 

IT 

we can parametrize the R-charges as 

R[Xi] = i?i2, R[X 2 ] = R41, R[X 3 ] = -R34, R[X 4 ] = R 23 , (3.27) 

in the electric theory and 



R[Y X ] = i?34, R[Y 2 ] = R23, R[Y 3 ] = R 12 , R[Y A ] = R A1 , 
R[Z X ] = R 42 , R[Z 2 } = R31, R[Z%\ = R-2A-, R[Za\ = R13 



(3.28) 



in the magnetic theory. Note that there are relations 

R[Xx] = R[Y 3 ], R[X 2 ] = R[Y 4 ], R[X 3 ] = R[Y X ], R[X 4 ] = R[Y 2 ] , (3.29) 
and 



R[Z,} = R[X ± ] + R[X 2 ], R[Z 2 ] = R[X 2 ] + R[X Z ] 
R[Z 3 ] = R[X 3 ] + R[X 4 ], R[Z 4 ] = R[X 4 ] + R[X X ] 



(3.30) 



The last four equations are natural since Z{ is the meson composed of two electric quarks 
Xt's. 
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Figure 11: The labeling scheme for the bifundamental fields and zig-zag paths for Seiberg 
dual theories. 



Let us compute the superconformal index. The standard Seiberg duality claims that 
electric and magnetic theories flow to the same IR fixed point. In the electric theory we 
have SU(N) gauge theory with Nf = 2N flavors. This theory has in general SU(Nf) x 
SU(Nf) flavor symmetries, but in our case we choose SU(N) 41 subgroup and denote the 
chemical potentials by (si,t, s 2j i) f° r fundamental flavors and (ti^, t 2ti ) for antifundamental 
flavors (i = 1, . . . , N), satisfying Yii s k,i = IL^M = 1 for & = 1, 2. From the rule of the 
previous subsection the index is 



Ie(p, q) 



K(p, q) 



N-l 



N\ 



TT dz j EL^i Ul=i T (SaZj-,p, q)T{T a z i 1 ; p, q) 



27rizj 

i=l J 



Y[ i ^ j T{z i z j ;p,q) 



(3.3i; 



for 



R[x 3 ] 



{Sa} = {(pq) 2 ' s lti u 1 u 2 1 , (pq) 2 3 s 2 , i u 3 u 4 1 } , 

R[X 2 ] _-. i fl[X 4 ] _ 1 _ 1 



(3.32) 



The magnetic theory is again has SU(N) gauge group with 2N flavors, and the index 
contains contributions from mesons: 



re(p, q) 



N-l 



2N 



l M 



N\ 



[r(u a , b ;p,q) / ; 



a,b=l 



N-l 



dz J IL=i IL=i T ( S aZj]P, q)T(T aZi 1 ;p,q) 



1=1 J 



Y\ i ^i.z i z j l ]p,q) 



for 



i S a} = {{pq) 2 s hi u 3 u 4 , (pq) 2 s %i u x u 2 } 



R[Y 4 ] 



{T a \ = {(pq) 2 t lti u 2 u 3 , (pq) 2 t 2 ,iU i u 1 } . 



(3.33) 



(3.34) 
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and 



R[gl] _-, _-, -R[Z 2 ] i _i 

R[Z 3 ] _- _- _R[Z 4 ] _- _- 

(pg)^~s7i*2,i«2M4 > (P9) 2 s l,^2,jWiU 3 } . 



(3.35) 



What we want to prove is the equivalence of the two expressions 

lE(p,q) = hi(p,q) • (3-36) 



As in |J8], we will establish this with the help of a remarkable identity of elliptic hy- 
pergeometric functions proven by [19], which is one of the many identifies studied, for 
example, in [50| IBTj. 

In order to apply [39], we need to check the "balancing condition". To state this, 
define 



H[*ll+fl[X 3 ] 



-1„. „.-l 



S: = (Y[S a ) 1 ' N =(pq)~ 



(3.37) 



a 

The balancing condition states that 

ST = pq , (3.38) 

which follows from the expression for S and T above. In this case we have Ig = J^-, 
provided 

S a = S/S a , f a = T/T a , U a , b = S a T b . (3.39) 

These conditions follow from the equalities ( I3.29p . (13. 30 p . 

It should also be kept in mind that to establish the identify above ( 13. 36ft we do not 
need to know the exact values of the R-charges determined from a-maximization; the 
equality of the index holds before a-maximization. 

3.3 Z-invariant Spin System 

We next show that the superconformal index for our quiver gauge theory can be refor- 
mulated as a classical spin model on a lattice in T 2 . 

The basic idea is simple. Let us first regard the integral variables a v = (a\, . . . , a v N ) 
as a ^-component continuous spin variable at vertex v G V (recall that these parameters 
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are related to z v by z v = e l(Jv ). These variables are circle valued with period 2ir, and in 
the SU(N) case the N components satisfy a constraint J2i=i o\ = 0. 

Next we need the Boltzmann weights associated with spin configurations. This is 
determined from the 1-loop determinant 7^ ect and I^nai- 



4d = Z spin = lUdv e ) e -E e6B ^H-E,^^k] 5 (3.40) 




where we introduced a new expression 

e = -'vector' e e = -'chiral ■ (3-41) 

In this language, I" ect is regarded as the self-intersection of the spins at position v G 
V, and the latter, 1^^, is the nearest-neighbor interaction of the spins at positions 
s(e),t(e) G V. 

To some readers this might look like a trivial re-naming of what we already know. 
However, first note that it is not true for general quiver gauge theories that resulting spin 
system is realized on T 2 ; the assumption of toric Calabi-Yau dual was crucial for this fact. 
Second, what is surprising about this spin system is that it is integrable. One simple 
way to see this is that the invariance of the index under the move II is the invariance of 
the partition function of the spin system under the double Yang-Baxter move, and the 
invariance under the Yang-Baxter move is one form of the integrability of the model. The 
chemical potentials of the 4d index (p, q, u) are regarded as the rapidity variables. In fact, 
the spin system constructed above coincides with the spin system studied by Bazhanov 
and Sergeev (hereafter BS) [Til [T2] . modulo some differences mentioned below. 

One technical difference is that the BS model are defined on the plane M 2 , whereas 
our model is defined on T 2 . Another more essential difference is that the integrable 
models have an invariance under a single Yang-Baxter move, whereas our index has an 
invariance only under the double Yang-Baxter move. 

Let us here explain this difference in more detail. In the BS model, to realize a single 
Yang-Baxter move we have to abandon the admissibility condition on zig-zag paths. Even 
in this case we could still choose a checkerboard pattern for the faces, corresponding to 
the colored and uncolored faces in our previous discussion. More concretely, in the brane 
realization explained in section 12.31 regions with (N, fc)-brane with k even (k odd) are 
colored (uncolored). We associate the spin variables a V)i to uncolored regions, where 
i = 1, . . . ,N and v is the label for the uncolored region. 

To define the Boltzmann weight, first we define the self energy e~ £v ^ as before. The 
definition of the nearest neighbor interaction e~ £e ^ is more tricky, since in this general 
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case we have two different types of edges, shown in Figure [12j For the two possibilities 
BS model assigns two different weights, which was denoted by W and W. 

What BS has shown is that for a judicious choice of £® s and Sf s , the model is 
integrable and is invariant under moves shown in Figure El including the single Yang- 
Baxter move. For N = 2 the weight satisfies the star-triangle relation, and this reduces 
to Spiridonov's formula for the elliptic hypergeometric function [52]. For N > 2, the 
start-star relation is still a conjecture, although there is non-trivial evidence from power 
series expansion, see [TTj . 





(a) (b) 

Figure 12: In admissible configuration of zig-zag paths we only have (a) as a possibility. 
However, if we lift admissibility condition, there are two possible types of intersections, 
and correspondingly we need two different weights. Red lines represents zig-zag paths, 
and gray region represents colored region. 

How is this related to our index? First, let us specialize the BS model to the case 
with admissible zig-zag paths. Then the case (b) in Figure [12] does not arise. The claim 
is that in this case the Boltzmann weights of BS model coincides with that of the 4d 
index, up to a spin-independent (but rapidity-dependent) overall normalization of the 
partition function. 

Let us see this explicitly for N = 2. In this case we can write z v = (z v , z~ l ) and the 
weights are 

^chiralO^e), Zf(e)) = F((pq) 1/2 ^j^^W^fe)^' «) ' ( 3 " 42 ) 

where here we meant the product of four terms with plus/minus signs. We also have (see 
(EO|) and (jMD^ 



c*(*) = (pvr 1/8 8i (vk 1/2 ) *i (vl ql/2 ) ■ (3 - 43) 
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Contrary to |llj we have included the factor 1/N\ in the measure, not in I v 
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These are identified with the weights function W Q and S of [llj 11 ! . except by an overall 
spin- independent normalization factor for 7 c hi ra i- 

The subtle difference in the normalization arises since statistical physicists and gauge 
theorists have different motivations. First, BS model requires a stronger invariance, i.e., 
invariance under moves violating admissibility conditions (Figure |SJ), and the specific 
normalization was chosen to serve that purpose. This is not necessary for 4d index. 
Second, 4d index has a very specific normalization build into the definition, whereas as 
a statistical mechanical problem we can multiply /vector by a constant depending on the 
chemical potentials p, q and u. For example, the two moves of Figure |6] preserves the 
number of the vertices of the quiver, and any overall multiplicative factor associated 
with the vertex hence preserves the invariance under the two moves. In fact, contrary to 
|llj . we are not normalizing the values of the partition function in the thermodynamic 
limit. 



As we explained in section I2.3[ the reason we stick to the admissibility condition is 
that we do not have a known Lagrangian for brane configurations for general D5/NS5 
brane configurations. The fact that we can extend the integrable spin system to non- 
admissible configuration could be suggesting that the Seiberg duality extends to theories 
on multiplet NS5-branes. 

Comments 

Several comments are now in order. 

First, the statistical mechanical model here is invariant under the moves in Figure El 
and has Z-invariance [531 G3] (studied originally by Baxter in the context of eight- vertex 
models). 

Second, in general the Boltzmann weight for the model is not real and positive [T2] . 
which might look problematic as a statistical mechanical model. However, in our context 
we are computing the index and the answer necessarily comes with negative contributions. 

Third, note that the resulting spin model is in general chiral, except for the case 
N = 2, where there is no distinction between fundamental and anti-fundamental rep- 
resentation. This explains why the mult i- component spin models of [12] requires the 
chirality (and therefore extra combinatorial data), as opposed to the N = 2 case of [TT] . 

Fourth, in the integrable models there is a duality transformation exchanging the 
graph Q with its dual Q*, which is a generalization of duality of the Using model ex- 
changing high-temperature expansion with the low-temperature expansion. In our gauge 



L The parameters p, q here are denoted by p 2 , q 2 in jllj . 
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theory context this duality is broken; for example Q* is a bipartite graph, but Q is not 
in general. 

Fifth, in the analysis to this point we have concentrated on the 4d index on S 1 x S 3 . 
However, as is clear from the argument above the specific choice of this observable in itself 
does not matter — other quantities will serve the same purpose as long as (1) they are 
invariant under the Seiberg duality and integrating out massive matters (2) the quantity 
takes a factorized form as in ( 13.101) . 

As an example of such quantity, we could use the 4d lens space index computed in [29] . 
This is a index on S* 1 x L(p, q), where L(p, q), with p and q coprime integers, is the lens 
space (Z p orbifold of S 3 ). The complication for this case is that we have a set of integers, 
which parametrize the discrete Wilson line. This means that the corresponding spin chain 
has integer as well as continuous labels, and could not be interpreted as a conventional 
statistical mechanical model. This lens space index has a reduction to 3d superconformal 
index in the p — > oo limit [29], which suggests a generalization of the contents of the next 
section concerning 3d S 3 partition function to the 3d superconformal index. 

Finally, let us make a brief comment on the comparison with the existing literature. 
The invariance of the index under the Seiberg duality for quiver gauge theories has been 
studied for some examples in [55]. The results of this section apply to more general quiver 
gauge theories dual to an arbitrary toric Calabi-Yau manifold. Note also the relation of 
the 4d index and the statistical mechanical model is known to experts, for example in 
[52] . However, to the best of my knowledge the correspondence with the supersymmetric 
quiver gauge theories has never been worked out in detail, and as we have seen there are 
in fact some subtle differences between the statistical mechanical models in the literature 
and the spin system defined from the 4d superconformal index. Moreover, as we will see 
in the next section we will see that this insight leads to a remarkable connection between 
3d M = 2 gauge theories and the geometry of 3d hyperbolic 3-manifolds. 

4 Reduction to 3d and 2d 

4.1 Reduction to 3d 

The spin system defined in the previous section is based on a rather general solution 
of the Yang-Baxter-type equation, and several known integrable models, including the 
Kashiwara-Miwa model [56] and the chiral Potts model [571 [58j [59] (and their sl^ gener- 
alizations), arise as a specialization/limit of the solution. Therefore it is natural to ask 
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whether there are some natural limits which is of direct interest in the context of gauge 
theory. 

In gauge theory there is one obvious limit. Since we have a 4d theory on S 1 x S 3 , we 
could dimensionally reduce along the S 1 . The resulting theory has 3d M = 2 symmetry, 
and by flowing to the IR we obtain a new 3d SCFT. We always need to keep in mind 
that there is non-trivial RG flow involved in this process; in particular, the anomalous 
dimensions in 4d and in 3d are in general different. This reduction is also a natural in 
terms of integrable models — it is simply the high-temperature limit of the theory. 

The effect of this S l reduction on the superconformal index has been studied in 
[60| |6T| EHJ [62] , and we obtain a 3d partition function on ellipsoid S$ [63] (generalizing 
the earlier results for S% =1 [&3J [65], 166] ) 

Sl = {(z u z 2 ) e C 2 | 6 2 N 2 + r 2 |z 2 | 2 = 1} . (4.1) 

To see this explicitly, first note that when the thermal S 1 shrinks all the chemical 
potentials p, q, u go to 1, but we could keep the ratio fixed and finiteo 

p = e -^+v) i q = e -W-v) i Ui = e- p,H . (4.3) 

For a 1-loop determinant of the chiralmultiplet we have (see ( 13. 18[) ) 

I _ e -/3[-*(o" s ( e ),j-o- i{e)ifc )-(fl e -l+/i i(e) -Ati i ( e ))+(-m+i)(l+»7)+(n+|)(l-r ; )] 

chiral-J.1 11 " - /3 r i((Ts(e)j _ (Jt(e)ifc)+{Re _ 1+AtL(e) _ Mij(e))+(m+ l )(1+77)+(n+ l )(1 _ ??) l -l • I 
j,k m,n>0 L c 

In the limit /3 — > we regularize the expression 1 — e~^ x to be 

1 — e~ /3x 
[x]/3 ■= — p ^ x ' as /3 ^ , (4.5) 

and Ichirai reduce to Z^ hival , which is given by 

Z e "TT TT li^f^hi ~ a *( e )' fc ) ~ ( re - 1) + (™+ j)(l + v) + (n + |)(! ~ V) 

Chiral iiiiio ^(^ - a ^ + ^ -V + (m+ |)(1 + V) + (n + |)(1 " V) 

TT TT j ("^(^(ejj ~ t7f(e),fe) ~ Q]e ~ jj) + ( m + j)& + (^ + |)& -1 

_ /ii m Jo f H^WJ " *t(e),*) + ^ - !)) + (^ + D 6 + (" + D 6 ^ 
N / Q \ 

= JJ s fe f <7 s(e))i - d- t(e))fe + L_ (1 - r e ) J , (4.6) 
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In t, y, u variables this is to take 

t = e-< 3 / 3 , 77 = e - "", m = e~ te , (4.2) 

which coincides with the limit taken in [HI] for 4d A/" = 2 theories. 
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where we defined 



(4.7) 



1 — T) 

and Sb(z), the quantum dilogarithm function defined in appendix, is used to regularize 
the infinite product (1A.20D . In the last line we re-defined a by a factor of ^: a := ^a. 
We also defined the R-charge r e by 

r e = R e + /x L(e) - fi R{e) . (4.8) 

Here the parameters fa play the role of the real mass parameters obtained by weakly 
gauging the i-th global symmetry. The equation f!4.8[) shows that the real mass param- 
eters has the effect of changing the anomalous dimension (cf. [67]). This is consistent 
with the Romelsberger's prescription explained in section 13.11 

Similar analysis shows (after we regularize the divergence coming from K(p,q)) 

J vcct -> Z vect = II S M ^ ~ ^ + ~2 ) Sb ( ~(° v ' i ~ ° v >i} + ~~2 ) ' 

= J J 4 sinh 6 (a^^ - a^) sinhfr -1 (& v>i - & Vij ) . (4.9) 

i<j 

We will drop the overall multiplicative constant 4 in the following. Note that the invari- 
ance of the index under y -h- y^ 1 (j3.2p is translated into the invariance b -H- 6 _1 . The full 
partition is still written as a matrix integral of these 1-loop determinants 

Z 3d = J [j[daA (I] Sector ) ( I! ^chiral ) , (4-10) 

which is precisely the 3d partition function on an ellipsoid. Note that in our story not all 
the real mass parameters of the theory are turned on. This is because we started from 
the 4d theory; some 3d global symmetries are anomalous in 4d and cannot be included 
in the 4d index. 

Our 4d superconformal index has an invariance with respect to the two fundamental 
moves. Our 3d partition function to have the same property, since the latter is simply 
the limit of the former. For example, in the same notation as in f)3.25p . the effect of move 
I is represented by 

JJ s b i a V j - 6v,j + — (1 - r e ) J JJ s b l <v,j - &vj + -y (1 - ( 2 - r e )) J . (4.11) 

i,j ^ ' i,j 

This is trivial due to the relation (1A.22J) . Of course, this is simply the limit of the relation 



(JMSD. 
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4.2 Higgsing 

To make contact with the geometry of hyperbolic 3-manifolds, we need to do one more 
reduction. For clarify let us specialize to the case N = 2 from here on. 

As stated in the introduction, we give a VEV to the vector multiplet scalar for the 
diagonal gauge group t/(l) diag 

[/(lUag C SU(2)W ; (4.12) 

and send the VEV to infinity. 

More concretely, let us write the vector multiplet scalar for the £/(l)diag symmetry as 
<j, and Higgs the £/(l)dia g gauge symmetry by sending a to infinity: 

i_old . ".new , ~ \ v »new r\ Z- > ~„ (A 19\ 

<j v — > a v + a, > a v = (J, a — > oo . (.4.13) 

vev 

In this limit, the vector multiplet 1-loop determinant diverges to a constant independent 
of a v , and after regularization we have 

^vector ~~ ^ -" vector = ^ ■ (.4.14) 

This is simply the 1-loop determinant for the Abelian gauge theory, which means that 
the gauge group after Higgsing is simply given by 



\U(1\ 



tf(l)dh« • (4-15) 



Lvev 
For the chiralmultiplet, the 1-loop determinant in (14. 6 j) for N = 2 reads 



Sfe (cT s(e ) + a t (e) ~ Z -yK) S b (cr s(e ) - & t (e) ~ ^rf) 



Z chiral - — — ; iO~Z\ ~ ^ iQ~X ' [4.LQ) 



where we defined 

r* = 1 - r e . (4.17) 

In the limit (14.131) , two out of the four factors stays the same 

s b (±(a s ( e ) - dt(6)) + yr*J ->■ s b (±(a s{e) - a t{e) ) + l —r* e J , (4.18) 
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whereas the other two factors behava 13 ! 

s b (±(a s{e) + a t{e) ) + ^|r e *) -+ e ±¥(±(^(e)+^))+^e) 2 . (4.19) 

After subtracting the divergences, most of the finite part cancels due to the ± sign, and 
the only remaining term is 



However, this cancels out when we sum over e G E, because 

e£_E v£E e: around v 

where we used (14.131) and that (thanks to (I2.18P ) 

E < = 2 - ( 4 - 2 °) 

e: around v 

The chiralmultiplet 1-loop determinant therefore becomes 

7 _, 7 / _ s b (a s{e) - a t{e) + §r* e ) 

^chiral — > ^chiral ~ / iO I\ ' l 4 - /J J 

Sfe (o- s(e ) - a t{e) - -fr*) 

The total partition function still takes the form (I4.10p . except that the 1-loop deter- 
minants are replaced by Sector) ^chirai- Using (see (IA.17|) ) 



i2, , 

e b (s - l -fr) 



s b {s +^r)_ = cQsnr e b {s + ^ r 



we see that this coincides with the Boltzmann weight Wq(s) of the Faddeev-Volkov model 
[TBI [T4l [15] in [17] under the parameter identification between their rapidity 9* e lj and 
the R-charge r* (recall (12. 14ft ) 

e: = -Krl e e :=7i-e* e = 7rr e , (4.23) 

up to an overall normalization 15 !. Again, the overall normalization worked out in [T7] is 
not necessary for our purpose. This two-step reduction from the solution of the star- 
triangle relation in [UJ to the Faddeev-Volkov model has been explained in a different 
manner in 1521. 



13 The appearance of the exponential of the quadratic term represents the parity anomaly. This cancels 
out in our final expression. 

14 This is denoted by 9„ in [T7j. 

15 Our function e& is denoted by tp in [17] . 
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The Faddeev-Volkov model is realizes Virasoro algebra on the lattice, see also the 
subsequent formulation of the discrete Liouville theory [68] . As explained in [18] , Liouville 
theory could be thought of as a boundary theory of SL(2, M.) Chern-Simons theory, and 
plays crucial roles in the relation with 3d M = 2 gauge theories. This suggests a more 
direct relation between the 2d spin system and the 3d hyperbolic geometry. In the story 
of [18] the geometric picture simplifies in the semi-classical limit of the Chern-Simons 
theory [20]. On the gauge theory side this limit is translated into the limit b — > 0, where 
S% in ( 14. ip degenerates to M. 2 x Si (the radius of Si is small, and is given by b). In the 
following we are going to take exactly the same limit. It is worth pointing out that the 
same limit of the Faddeev-Volkov model has been analyzed in the seminal paper [17] , 
which also studied the connection with circle patterns. 

4.3 Further Reduction to 2d 

Let us dimensionally reduce our theory further to 2d, by taking the b — > limit of the 
S% partition function. The resulting theory has 2d M = (2, 2) supersymmetry. After the 
dimensional reduction 2d the vector multiplet scalar a is complexified due to the Wilson 
lines along Si, and a 3d real mass parameter reduces to a twisted massed in 2d. 

In this limit, quantum dilogarithms reduce to classical dilogarithms (1A.24I) . and we 
have 



Z' 3d ->• Z 2d = / da exp 



" W 2d (a) 



7Tb 2 

where the effective twisted superpotential W2d is given by 



(4.24) 



W 2d ((T) = W 2d (p) = -- J2 [HPs(e) ~ Pt(e) + ie* e ) ~ I (p S (e) " Pt(e) ~ Z^)] • (4.25) 

e€E 

We have scaled the variable a by 2nb (this comes the radius of Si as in the standard 
dimensional reduction) 

p e = 2nba e =nQba e , (4.26) 



and l(x) is defined in ( 1A.25J) in Appendix. Because l(x + 2iri) = l(x), we can regard 9 e 



as defined modulo 2ttZ. This is consistent with the discussion in section 12.41 Note also 
the appearance of the classical dilogarithm represents the sum over the contributions of 
the KK modes along Si [69]R . 



16 There is a side remark about the integrable structure. In references |69l 170] . the saddle point 
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The classical vacuum of this theory is determined from the saddle point equation of 
the twisted superpotential 

exp {-^r) = l ■ (427) 

In our case, this is given by the "cross ratio equation" (which are closely related with the 
Hirota difference equations) 

1= TT cosh (p 8(e) - p t{e) + 16*) yr cosh (p 8(e) - p t{e) - i9* e ) 

4.4 Circle Patterns and Tetrahedra 

Surprisingly, the effective twisted superpotential (J4.25P as well as the saddle point equa- 
tion (14.281) has a beautiful geometrical reformulation. For this purpose we need to add 
extra degrees of freedom to the bipartite graph. In most of our discussion in previous 
sections, the bipartite graph only contains combinatorial data — we do not need to spec- 
ify the length of the edges (the only exception is the part where we discuss isoradial 
condition )o. Henceforth we parametrize the length and angles of the bipartite graph 
and zig-zag paths in the following manner. 

Let us assign a parameter p v at each vertex v G V, and draw a circle with radius 
r v := e Pv centered at v. The quantity p v will be identified with the vector multiplet 
scalar ( 14 . 2 6 [) momentarily, and hence we use the same symbol. We moreover impose the 
condition that when v and v' are adjacent at an edge e the two circles intersect at an 
angle 9*. This angle will again be identified #* introduced previously in ( 14.23ft . Let tp v ^ e 
to be an angle around vertex v (Figure [TBI , and elementary trigonometry shows that 

= e* + e^i = „. cosh(p v - Pv , + z6;) 

e p v + e p v/ -id* e C0sh(p„ - /V -i0* e ) ' 



equation of effective twisted superpotential of a 2d theory is identified with the Bethe Ansatz equation 
of an associated integrable model. It is interesting to imagine if this integrable structure (which in turn 
is conjectured to be equivalent to the integrable structure of the dimer model [26]) could be related to 
the integrable structure of the spin system mentioned above. However, it should be kept in mind that at 
least the two should not simply be identified directly; the two spin systems are rather different. Another 
difference is that the rapidities in our spin system are the chemical potentials for our 4d index, whereas 
those in the context of [69j [70] are identified with the vector multiplet scalars. 

1 In the context of brane configurations studied in section [2~3l thcse quantitative datum do matter. For 
example, the area of the region of the (N, 0)-brane represents the inverse square of the gauge coupling 
constant. However, the author is not aware of direct brane interpretation of the geometric construction 
presented in this section. 
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see Figure [T3j We also have 



<Pv,e + fv',e = 2#* 



(4.30) 



Note that the right hand side of this equation takes the same form as the expression in 
(14.28)) . Thus the cross-ratio equation is the statement that the angles around v sum up 

to 2tt 



2 4 We = 27T , 
e: around v 

where we used (see ffl~20"]) and (H~23j) ) 



J2 K = ° (mod27rZ) 



(4.31) 



(4.32) 



e: around v 



vertices/of the bipartite graph 




vertices of the quiver diagram 



Figure 13: The two circles of radius r v and r v >, each centered at vertex v and v' connected 
with an edge e, intersect at an angle 9*. This determines the angles tp V!e ,(p v ' !e . The red 
quadrilateral is sometimes called a kite, and is a rhombus for an isoradial circle pattern. 

Summarizing, the saddle point equation (14.281) is reformulated as a geometric condi- 
tion on a set of circles intersecting at the vertices of the bipartite graph. In the literature 
a set of such circles is called a circle pattern 18 !, and the deformation of a circle pattern 
represents the discrete analog of conformal transformation. 



IS 



In the case that 6* = ? for all e E E, a circle pattern reduces to a circle packing together with its 



dual. 
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Figure 14: Our polyhedron M is defined as a complement of all the hemispheres. 

For us, the importance of this geometric reformulation is that it gives the 3-manifold 
we are aftecj. To explain this, let us regard T 2 as part of the boundary M 2 of H 3 , and for 
each circle consider a hemisphere centered at v and intersecting the boundary M? with 
the circle around v. Let us denote the region above all the hemispheres by M (Figure 
IT4"|) . Note that we still keep the torus identification on the boundary. 

This could be equivalently described as follows. For each face v G F* = V of the 
bipartite graph Q* you can associate an ideal polyhedron A„ whose vertices are the 
vertices of Q* at the face and at infinity. Then the polygon is defined as the union of all 
these ideal polyhedra, with neighboring faces glued together: 



M 



IK- 



(4.33) 



vev 



This is a hyperbolic 3-manifold (in the standard metric induced from that of H 3 ), and 
has a geodesic boundary. In particular its hyperbolic volume is computed to be the sum 
of the volume of A„, and depends non-trivially on p v and 9*. 

The manifold M can be decomposed into hexahedra (Figure [T5]) . which in turn could 
be divided into two non-ideal tetrahedra. As is clear from this definition, the gluing 
condition of the hexahedra at the vertices of the quiver diagram is exactly the condition 
(14.3 ip . ensuring that the total angle around an edge of the polyhedron is 2n. We have 
therefore shown that the vacuum equations of the 2d gauge theory coincides with the 
gluing conditions of the 3-manifold M. As discussed in [1] , the existence of the solution 
of the gluing condition is guaranteed by the 4d conditions on the R-charge ( 12. 151) . ( 12. 17ft 
and the result of US]. 



19 Historically the fact that a circle pattern has to do with the geometry of 3-manifolds is known since 
long ago, see for example [7TJ Chapter 13]. 
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We can also directly compare the value of the twisted superpotential with the volume 
of the hyperbolic 3-manifold. We show that 

W^)|_ - WmL, = Vol[M](a)| glumg - Vol[M ] , (4.34) 

where the quantity Vol[M ], to be defined below, is independent of p v , and we have 
evaluated the both sides at the values of p v determined from the vacuum equations or 
gluing conditions. To evaluate the volume, we use the decomposition of M into hexahedra 
(Figure [T5|) . The volume of the hexahedron at an edge e is given by [721 Lemma 4.2] 

lt-» = " (^p) + * {^T) (4-35) 



l(Ps(e) - Pt(e) +i0* e ) -l(i6* e ) 

i(p S (e)-Pt(e)-ie* e )-i(-ie:) 



+ 2L 

+ 7(Vs(e),e - <Pt(e),e)(/>s(e) - Pt(e)) , 

where in the last line we used (1A.14I) . ( 1A.27I) in Appendix. The total volume of M is 
then given by summing over the edge e: 

™[M] | gluing = E Kexahed J glumg = V^| rfntag " W 2d | ct=0 + Vol[M ] , (4.37) 

where we used that fact that the last term in (I4.36[) cancels out at the saddle point (see 
( 1413J1 . Km . (0L32D), and we defined 



Vol[M ]:=E 2jI (|) • ( 438 ) 



As we will see shortly this is the value of Vol[M] when the circle pattern is isoradial. 
This proves (14.34)) . 

Isoradiality 

In general, it is a difficult problem to analytically solve the cross ratio equations 
(I4.28p . Considerable simplification occurs when the graph is isoradial. Isoradial condition 
(section I2.2J) states that it is geometrically possible to choose p v such that 

p v = p independent of v . (4.39) 

This is particularly useful for us because in this case (I4.29P simplifies to 

ip e = 9* e (mod 2ttZ) , (4.40) 
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hexahedron 




Figure 15: Our polygon M could be decomposed into hexahedra, each of which is asso- 
ciated with an edge of the bipartite graph. The projection of this hexahedron gives the 
red quadrilateral in Figure | 



which automatically solves (14.311) (see (14.321) ). There is an enhanced S^d^l) global 
symmetry at this vacuum. 

In this vacuum the volume of M reduces to (14.381) . and the critical value of the twisted 
superpotential vanishes, 

W 2d (a) = . (4.41) 

isoradial 

Note that the isoradiality is not preserved under the two fundamental moves of Figure 
EJ However, the critical value of the twisted superpotential is still preserved under the 
fundamental moves. 



Hyperbolic Circle Patterns 

To this point we started with a geometric realization of the dimensionally reduced 
theory after the Higgsing. This raises a natural question: what if we start with the 3d 
theory before Higgsing and dimensionally reduce to 2d? Does this theory has geometric 
reformulation? The answer to the latter question is affirmative. Let us explain this 
briefly. 

Let us go back to the expression of the partition function in section 14.11 an d again let 
us take iV = 2, in which case the chiralmultiplet 1-loop determinant is given in ( I4.16p . 
After the dimensional reduction as in section 14731 (but without going through the Higgsing 
in section |4T2|) . the 2d vacuum equation is written as 



£ 



¥>«,. 



2vr, 



(4.42) 



e: around v 
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where tp VtE is defined by 



D 2iifi Vt , 



\ _|_ gf„', e -p»,e-tfe ^ I e P v > , e +Pv,e+iO e 



(4.43) 



where the edge e connects the vertex v with another v'. The equation (I4.43P is a trigono- 
metric relation for the angle placed on Figure [131 but now placed on a hyperbolic plane 
H 2 instead of R 2 [72]. From this viewpoint the Higgsing in section 14.21 is a process of 
replacing circle patterns in H 2 by those in M 2 . 



4.5 Comments on the Geometry of the 3-manifold 

In the previous subsection our 3-manifold is obtained by gluing ideal polyhedra. Here 
are preliminary remarks on the geometry of this 3-manifold; more details are currently 
under investigation. 

M5-brane Realizations 

Let us comment on the brane realization of our 3d M = 2 theory. We begin with the 



brane configuration in section I2.3[ and in order to dimensionally reduce the theory to 3d 
we compactify and T-dualize along the 3-direction. Then we have a D4/NS5 system in 
type IIA: 



(4.44) 
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This can be lifted to M-theory by including an extra 11-th direction: 
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(4.45) 



The resulting system is the intersection of two different types of M5-branes, where N M5- 
branes (coming from D5-branes in type IIB) wraps T 3 and another M5-brane (coming 
from the NS5-brane) wraps a special Lagrangian sub manifold M' in IR 3 x T 3 . Note 
that M' and T 3 , two different types of M5-branes, intersect along a one dimensional 
subspaca 20 !. This brane configuration tells us that M' is the Riemann surface fibered 
over the S 1 ; of course, this S 1 is part of the S 1 x S 3 on which we compactified 4d gauge 



3 Similar brane configuration has been studied in [75] , 
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theory. This places a strong constraint on the geometry of M' . It is a future problem to 
clarify the precise relation between M' and our 3-manifold M. 

M As a Link Complement 

Our 3-manifold M is a link complement, at least in some special situations. Let us 
illustrate this with a simple example, following [HJ Appendix] and [75] . 

Suppose that we start with a bipartite graph shown in Figure [16] (a). The corre- 
sponding Calabi-Yau geometry is the canonical bundle over P 1 x P 1 . Suppose that the 
R-charge of all the fields are 1/2. The circle pattern then is a circle packing with its dual 
(Figure [TBI (a), green circles). The bipartite graph Q* have four faces (gauge groups), 
and we choose a checkerboard coloring of the faceo (Figure [TBI (b)). We first glue two 
colored faces. The resulting manifold has a boundary. We next prepare an identical copy 
of this partially glued polyhedron and then glue the two along the remaining uncolored 
faces. This gives a link complement, shown in Figure [TBI (d). 

More generally, this procedure works if (1) a circle pattern is given by a circle packing 
with its dual, i.e., when the R-charges are all canonical, (2) the vertices of the quiver 
diagram can be checkerboard colored and (3) choose a particular pairing of the black 
faces. For example, this works for all the orbifolds of the example in Figure [TB] with 
canonical R-charges. 

5 Dimers and BPS State Counting 

In this section we study the relations of our results to BPS state counting and topological 
string theory. 

5.1 Thermodynamic Limit of Dimers 

Let us first begin with the observation that the Legendre transform of the volume of 
the 3-manifold M when the circle pattern in isoradial coincides with the thermodynamic 
limit of a dimer partition function. 

Suppose that I have a dimer model on T 2 . We choose a parametrization of the weight 
e Ue for an edge e G E by an angle ("rhombus angle") 8*: 

e^ = 2sm^.. (5.1) 

2 v ; 



21 This coloring is different from the coloring described in section l!01 
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identify 

the two faces 
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nt n 



(d) 



each identified with the boundary 
of another copy of M 



Figure 16: From two identical copies of our 3-manifold M we obtain the link complement 
shown in (d). We first start with a circle packing and its dual (a). We then choose a 
checkerboard coloring of the faces (b), and then pairwise glue the colored faces. This 
leaves a polyhedron with uncolored faces as boundaries, together with its copy. We finally 
glue the two polyhedra along the uncolored faces, and obtain a link complement in (d). 
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Note that this is the logarithm of the length of the edge e in Figure [13j The partition 
function of this dimer model is defined by 

m: perfect matching e€m 

where the sum is over perfect matchings, i.e., subsets of edges of the bipartite graph 
containing each vertex once and exactly once. 

Let us consider the thermodynamic limit of this partition function. As is standard 
in statistical mechanics, this limit is taken by enlarging the size of the torus by n times 
in the two directions of the torus. Let Z n be the partition function of this model. In 
the large n limit, the free energy scales as n 2 , and we could extract a finite quantity by 
taking the scaling limit. In the case that the dimer is isoradial, this is computed to be 
[761 177^ 

27rlim logZ n = y](0:iog2 S in 1 f + 2Jl( 1 f ) ] =:V(u) . (5.3) 

This quantity is originally introduced as a normalized determinant of the discrete Dirac 
operator [76]. This is consistent with the fact that dimers are described by free fermions, 
and the dimer partition function is the determinant of the Dirac operator. 

The expression (J5.3P is not exactly the volume itself (I4.38|) . but its Legendre trans- 
form. In fact, there is a general analysis of the thermodynamic limit of the dimer model 
in [75], which is applicable to arbitrary (not necessarily isoradial) dimers. There it has 
been shown that the Legendre transformation of the thermodynamic limit of the dimer 
partition function (the left hand side of (|5.3J) ) is the Ronkin function R (explained in 
Appendix [Bj of the spectral curve of the dimer model, and is given by the Legendre 
transformation of the surface tension a: 

V» = ^> e + a(r) , (5.4) 

eG-B 

where 9* is the conjugate variable to u e . By definition 9* can be computed as the 
probability of the edge e to be chosen in a perfect matching, and is shown to be related 



to 9* by (EL]) m- 

This fact is powerful enough to determine cr(9*). In fact, from the property of the 
Legendre transformation we have 

^ = -„ = -l„ g 2sin§, (5.5) 



J Our 0* here differs from in [H E7J by a factor of 2. 
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and by integrating with respect to 8* (see ( IA.9J) ) we obtain 



a(r) = ^2Jl(|)=Vol[Af ]. (5.6) 

By combining (j5.4p and (15. 6p we obtain (j5.3p . In Appendix [B] we demonstrate (I5.3P by 
an explicit computation for a simple example. 

Summarizing, we have seen that the Legendre transform of the volume of the polygon 
M coincides with the thermodynamic limit of the dimer partition function. 

Comments on F-extremization 

Before proceeding to a next subsection there is a side remark on IR anomalous di- 
mensions. The combination ( 15. 3 j) can be rewritten as 

'9* 



0*log(2sin- +2JI 

V 7 • / x ■ n (5-7) 

,log(l-e^) + M^ 2 + iLi 2 (e 2 -) -^ 



-2tt 



where 6** = 27rz and the function inside the bracket is the function Z(z) defined in [65] , 
which appears in the 1-loop determinant for the S% =1 partition function. This means that 
the extremization of the dimer partition function coincides with the extremization of the 
integrand of the S^ =1 partition function in [65] , and hence the IR anomalous dimension 
in the saddle point approximation of the matrix modeo. The author is not aware of an 
immediate application of this fact, nevertheless it would be interesting to explore if there 
is any further implication of this observation. 

5.2 Relation with Topological Strings 

We have seen that our partition function coincides with the thermodynamic limit of the 
dimer partition function. The same dimer partition function appears in the physics of 
BPS state counting and topological string theory (see [79j |80] for summary). 

Let us consider type IIA string theory on the toric Calabi-Yau X&. We consider BPS 
bound states of D0/D2-branes wrapping 0/2-cycles of X& bound to a single D6-brane 
filling the whole X aI 24 I . Equivalently these bound states are 1/2 BPS particles of the 4d 
Af = 2 supersymmetric gauge theory. 



23 A caveat in this statement is that there is a 3d global symmetry which is anomalous and not 
included in the superconformal index in 4d. 

24 We can define the crystal melting model with D4-branes wrapping 4-cycles included, however its 
thermodynamic limit has some additional subtlety due to the existence of the "gas phase" [75] . 
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We define the BPS partition function as the generating function of the BPS degen- 
eracies of these bound state of D-branes (mathematically this is the generating function 
for the generalized Donaldson-Thomas invariants): 

ZEPs(g.,t) = J2 n(p ,PA)e- 9 ^ p °- tAPA , (5.8) 

P0,PA 

where po and pa are DO and D2-brane charges, respectively, and Q is the BPS degeneracy. 

In the case the Calabi-Yau 3- fold is toric, then the BPS partition function is known 
to be described by the partition function of the crystal melting model, or equivalently a 
dimer model on M? [81 ] |38 | 182] . The resulting partition function takes an infinite product 
form, and is a reduction of a square of the topological string partition function |83) 25 l . 

The thermodynamic limit of the dimer partition function in the previous subsection 
appears in the thermodynamic limit of this dimer model, i.e., in the limit g top — > 0. This 
limit has been analyzed in [75] . and the limit (after suitable regularization) is an integral 
of the Ronkin function R(x, y) of the spectral curve of the dimer model, which as shown 
in [31] coincides with the Riemann surface appearing in the mirror of the Calabi-Yau 
manifold (I2.10p . This leads to the identification 



'bps -*• exp 



v 

—j-/ o 

9s 



^to P ,o = / dxdy R(x, y) , (5.9) 



where J-to P ,o is the genus prepotential of the topological B-model on the mirror Xa, 
or equivalently of the topological A- model on X&, and the integral is over the amoeba 
defined in Appendix |B] Combining these with our result (I4.38P we see that an integral 
of the Legendre transform of the critical value of the effective twisted superpotential 
reproduces the genus topological string partition function, as stated in (II. 7\i . 



^to P ,o= fdxdyC\Vo\[M] , (5.10) 

J L isoradial- 



where £ represents the Legendre transformation described in section 15.11 

As a simple check of this relation, let us count the number of parameters. As we have 
seen, the twisted superpotential W2d(o"), and therefore Vol[M](cr), have d— 1 parameters. 
Because we integrate over two of them, we have d — 3 remaining parameters. This is the 
same as the number of compact P 1 in the geometry, and hence the number of Kahler 
moduli, for a toric Calabi-Yau 3-folds without compact 4-cycles. 



25 Intcrestingly, this partition function could also be expressed as a matrix model [841(853, which is again 
written in terms of dilogarithm functions and is similar to the matrix models considered in subsections 
14.11 and !4. 21 The spectral curve of this matrix model reproduces the mirror curve. 
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4d/ld correspondence 

We have seen that a reduction of the 4d superconformal index for quiver gauge theo- 
ries reproduces the topological string partition function on the dual geometry, the toric 
Calabi-Yau manifold. As we have explained above the two quantities are apparently of 
very different origin: one comes from 4d M = 1 superconformal quiver gauge theory and 
another from the counting of BPS particles in toric Calabi-Yau. This raises the obvious 
question whether there is a natural explanation for this correspondence. 

The simple answer is that the same bipartite graph (quiver and the superpotential) 
describes the two different physics, (1) 4d quiver gauge theory and (2) Id quantum 
mechanics on the 1/2 BPS particles inside 4d J\f = 2 theory 26 !. In this description, 
Seiberg duality in 4d mutates the quiver, which on the Id side is interpreted as the 
crossing of the wall of marginal stability in the moduli space [861 EZ]- Because the 4d 
index is invariant under Seiberg duality, we should recover a quantity which is not affected 
by wall crossing, and topological string partition function precisely satisfies this criterion. 

This correspondence between 4d and Id is strongly reminiscent of the correspondence 
[22] between the 3d M = 2 gauge theories realized as domain walls in 4d M = 2 theories 
and the 1/2 BPS particles in the same 4d M = 2 theory, where the counterpart of a 4d 
Seiberg duality is played by a 3d mirror symmetry. 

6 Concluding Remarks 

We conclude with a few open problems, besides those already mentioned in introduction. 
Our results, in particular the relation (II. ip . are reminiscent of the recently found 
relation between 4d M = 2 superconformal index and 2d TQFT [EH [89], ED]- There 
are many important differences between the two proposals; we have different class of 
4d SCFTs (for example, our theories have Lagrangians whereas their theories do not in 
general) with different amount of supersymmetry (A/ - = 1 versus M = 2). However, it 
is instructive to pursue the analogy between the two. For example, the counterpart of 
the decomposition of 4d M = 2 SCFT from the pants decomposition of the Riemann 
surface [91] is the decomposition of the 4d M = 1 SCFT from the partial resolution of 
the toric Calabi-Yau singularities, and the mirror of the resolution is indeed the pants 
decomposition of the mirror curve (12.81 ) . 



26 This was one of the key ingredients which lead to the construction of crystal melting model }82| . 
The subtle difference is that in (2) we have a D6-brane filling the whole Xa- This plays crucial roles in 
the wall crossing phenomena, but not for the consideration of this section. 
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The relation between the 2d spin system and the 3d SL(2) Chern-Simons theory 
deserves further study. Since 3d SL(2) Chern-Simons theory is closely related with 3d 
gravity, our 2d spin system could be regarded as a kind of spin network deconstructing 
3d gravity. 

In the discussion of 3-manifolds in this paper, the natural starting point was the 4d 
superconformal index. At the technical level this means that the quantum dilogarithm 
function arises as a limit of the elliptic gamma function. Since quantum dilogarithm 
function appears extensively in the geometry of 3-manifolds, wall crossing phenomena 
and 3d J\f = 2 gauge theories 27 !, it is natural to ask if we could lift all these to the level 
of elliptic gamma functions. For example, is there a pentagon relation for elliptic gamma 
functions? 

A somewhat related problem is to give a mathematical formulation of the supercon- 
formal index. Since our 4d theory is defined purely from the toric data there should 
be a mathematical formulation of the superconformal index in term of the geometry of 
the toric Calabi-Yau 3-fold, perhaps as an equivariant character of some moduli space 28 !. 
Pushing this further, there should be "categorified" versions of (II. II) . ( 11.61) . for example 

4d J\f = 1 SCFT on 5 3 "moduli space defined from X& ; \y. i.) 



for some Hilbert space % (similar categorified viewpoint has been advocated in [T8]). 
Our results in section [5] suggest a relation of this Hilbert space with the moduli space of 
ideal sheaves on the Calabi-Yau manifold. 
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A Special Functions 

In this appendix we summarize formulas for special function used in this paper. 
We define the elliptic gamma function T(z;p,q) bjo 

1 - z~ l q j+l p k+l 



r(*;p,g)= ' [ , / k • (A-i; 

- LJ - l — ZQ^V 



j,k=o H F 

This is the basic building block for the 4d superconformal index. This function satisfies 

T(z;p,q) = - q —^ r, (A.2) 

T(pq/z;p,q) 

and 

T(pz; p, q) = 6(z; q)T(z; p, q), T(qz; p, q) = 6(z; p)T(z; p, q) , (A.3) 

where 6(x;p) is the g-theta function defined by 

00; q) = (z; q)oo(q/z; q)^ , (A.4) 

with 

oo 

(a;q)oo = l[(l-aq k ) • (A.5) 

fc=0 

g-theta function is related to the Jacobi ^-function 

oo 

e^zlq) = 2q^ 4 simrz ]\(1 - q 2m )(l - 2cos27rzq 2m + q 4m ) , (A.6) 



m=l 



29 Beaware of the notational differences in comparison with the literature. For example, our notation 
here is different from that in pT| . 
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by 



0i(zk)=tg 1/ V;^)ooe-~e(e*~;^). 



(A.7) 



The equation ( 1A.3I) is a elliptic generalization of the familiar relation T(x + 1) = xr(x). 
From flA~2|) . flA~3|) . (COD it follows that 

1 1 



r 0(*;p)0(z;g). 



(A.* 



(1 - z )(l- z^)T(z;p, q)T{z- 1 ] p, q) (1 - zf 

As for the 3d and 2d partition functions, we need classical and quantum dilogarithms. 
The classical dilogarithm function (Euler dilogarithm) Li 2 (z) is defined by 



Li 2 (z) = - 



log(l - 1) 



dt = Y J 



?1=1 



n- 



(A.9) 



A related function, Lobachevsky function Jl(x), is defined by 

Jl(x) = - dt log 1 2 sin 1 1 . 

Jo 

This could be expressed in terms of the Euler dilogarithm to be 

2i Ji(0) = Li 2 {e 2ie ) - Li 2 (l) - 6(6 - vr) , 
where the famous formula by Euler states that 

Li 2 (l) = C(2) = y • 

By definition the Lobachevsky function is an odd function, Jl(— x) 
spondingly, Euler dilogarithm satisfies 

Li 2 (-e*) + Li 2 (-e-*) = -" 2 
which is consistent with (1A. 12[) . We can also show 



7T 
X' 

6 2 



JI 



IX 



JI - 
2 



- 2 m 



(A.10) 

(A.ll) 

(A.12) 
Jl(x). Corre- 

(A.13) 
(A.14) 



We also need the non-compact quantum dilogarithm functions s^z) and e^z) [T 
[M] . These functions are defined by 



s b {z) = exp 



1 f°° dw ( sin 2zw z 

i Jo w \2sinh.(bw) smh(w/b) w 



(A.15) 
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and 



e b (z) = exp 



oo+iO 



dw 



-i2zw 



(A.16) 



-oo+io w sinh(wb)smh(w/b) 

where the integration contour in ( 1A.16I) is chosen above the pole at w — 0. In both 
these expressions we require \Imz\ < |Imc&| for convergence at infinity. There is a simple 
relation between the two functions 



e b {z) 



e 2 e 



Sb{Z) 



(A.17) 



and we loosely refer to both functions as quantum dilogarithms 
The function s b (z) satisfies a difference equation 

ib ±l 



Sb z 



2 cosh {jib ±1 z^ s b ( z + 



ib ±l 



and hence 



iQ\ ( iQ 

z —J =4 sinh {jib~ 1 z) sinh (wbz) s b ( z + — 



(A.18) 



(A.19) 



We can use these equations to analytically continue s b (z),e b (z) to the whole complex 
plane. The position of poles and the zeros are represented as 



nmb + nb + tt — iz 
6 — 
^„ mb + nb 1 + % + iz 



(A.20) 



where the left hand side is to be interpreted as a regularization of the right hand side. 
By definition we have 



s b (z 



s\/ b {z) 



(A.21) 



and 



s b [z s b i-z 



The asymptotic limit is given by 



s b (z) ->■ e : 



sgn(z)- 



as z — > ±oo . 



(A.22) 



(A.23) 



The name "quantum dilogarithm" could be justified by the fact that in the classical 
limit b — > 0, the quantum dilogarithm reduces to the classical dilogarithm: 



e b {z) — > exp 



2tt6 2 



:-i)Li 2 (-e 



2nbz\ 



s b {z) —> exp 



l(2nbz) 



2ixb 2 



(A.24) 
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where 



l(z) 



1 TV 2 

U 2 (-e z ) + -z 2 + — 



[Li 2 (-e 2 ) - Li 2 (-e" 2 )] 



, 2 7T 

2JI — + - 
>2i 2 



This function is odd, and its derivative is given by 

l'(z) = i log 2 cosh 



(A.25) 



(A.26) 



In the main text we need the following equality 



( y) - * (y) + jPte. - ¥>i] 



where ipi, tp 2 are defined by 

1 + e^ 



3 *</>i 



1 + eP- 



-if) 



1 _|_ p-p+i8 

^ 2 = TT7^> ^ + ^ = 29. 



(A.27) 



(A.28) 



The equality (1A.27P can be proven with the help of the formula (this is a minor modifi- 
cation of [95], Proposition A]) 



U 2 (-e p+ie ) - U 2 (-e p - ie ) = 2% [Jl(0) - Jl(w) - Jl(9 - u) - cop] 



(A.29) 



with 



tan a; 



e p sin 9 
1 + eP cos 9 



i.e., e 



2/w 



1 + e p+ ^ 
1 + eP~ ie 



(A.30) 



B Hyperbolic Volume and Ronkin Function 

In this appendix we first recall why the thermodynamic limit of the dimer partition 
function is given by the Ronkin function of the characteristic polynomial of the dimer 
model. 

We then show in an explicit example that the Ronkin function is the Legendre trans- 
form of the hyperbolic volume of the polyhedron M. This equality is known in the 
literature, for example for C 3 [96] and for a certain parametrization for the canonical 
bundle over P 1 x P 1 [97]. 

Consider combinations of v e along closed non-trivial loops (a- and /3-cycles) in T 2 , 
and let us denote them by e x ,e y . These are the chemical potentials for the so-called 
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"height function" of the dimer model. Written in this variable, the partition function for 
the enlarged dimer Z n (defined in section 15. ip is simply given by 

Zn(e x ,ey)= n n Z l (z',w'), (B.l) 

z' n =e x w' n =ey 

where only the dependence of Z n with respect to x, y are explicitly shown here. By taking 
a logarithm and dividing by n 2 , the product becomes a Riemann sum and we have 

lim ±logZ n (e x ,e y ) = / ^-^log |Z x (eV, e y w')\ = R Zl (x,y) , (B.2) 

where for a Laurent polynomial P(z, w) we defined its Ronkin function Rp to be 

Rp(x,v) = 7^ / — — ]og\P(^z^w)\ . (B.3) 

(27T2) 2 J\ zHwl=1 Z W 

This is a convex function, and is linear outside the amoeba Ap, defined by 

A P = {{x,y) E R 2 \ 3(6, <P) e (M/2ttZ) 2 , P{e x+l \ e y+i *) = 0} . (B.4) 

In general it is not easy to work out the exact analytic expression for the Ronkin function. 
However, in many cases it is much simpler to compute its derivative by taking the residues 
of the integral. 

Let us check the equality of (15. 3D and (1B.2J) (times n) when A is the toric diagram 
for C 3 . The dimer has three edges, whose weights we denote by x, y, z. Correspondingly 
we have three rhombus angles a, f3, 7, which is related to x, y, z by 

e x = 2sina, e 2/ = 2sin/3, e z = 2sin7. 

and satisfying 

a + (3 + 7 = ix . 

These weights satisfy the triangular inequality (we take a, (3, 7 > 0) 

e x + e y > e z , e y + e z > e x , e z + e x > e y . (B.5) 

The dimer partition function is 

Z x = e x + e y + e z , (B.6) 

whose amoeba coincides with the region determined by (IB.5[) . We now show 



irR Zl (x, y, z) = ax + (3y + -fz + Ji(a) + Jl(/3) + Jlfr) . (B.7) 
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inside amoeba, where we defined the 3-variable version of Ronkin function by 

R p ( x ,y, Z ) = -±-<f ^^l 0g |p( e * U)eVe% )| . (B.8) 

(2«) 3 J luHvHwl=l u v w ' ' 

This coincides with the 2- variable Ronkin function -R 1+ j> + £ up to a constant factor. 

xx 

To check (IB .71) let us compute the derivative of the both sides of (IB. 71) . The derivative 
of the left hand side inside amoeba can be evaluated by taking residues [98] 

/ e 2(x-y) _ e 2(z-y) _ ]\ 

(x - y, z - y) = tv - cos x ( ^^ J (B.9) 



ORz 1 , s ORl +e x +e y 

71 a (^. V-. z ) = n - 
ox ox 


and hence by the law of 


cosines 




dR Zl 



it - 1 = ix — cos l (— cos a) = a . (B.10) 

ox 

The computation is similar for the y, ^-derivatives, and the result coincides with the 

x, y, ^-derivatives of the right hand side of (1B.7J) . Note that the amoeba coincides with 

the region determined by (IB. 51) . 
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